Deceleration of continuous molecular beams 
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A method for decelerating a continuous beam of neutral polar molecules is theoretically demon- 
strated. This method utilizes non-uniform, static electric fields and regions of adiabatic population 
transfer to generate a mechanical force that opposes the molecular beam's velocity. By coupling 
this technique with dissipative trap-loading, molecular densities > 10 11 cm~ 3 are possible. When 
used in combination with forced evaporative cooling the proposed method may represent a viable 
route to quantum degeneracy for a wide-class of molecular species. 

PACS numbers: 37.10.Mn, 37.20.+j 



The electric dipole-dipole interaction between polar 
molecules is fundamentally different from most interac- 
tions between ultracold atoms. While atomic interactions 
are typically isotropic and comparatively short-ranged, 
the dipolar interaction is strong, long-range, tunable and 
anisotropic. This interaction can lead to many novel and 
exciting phenomena, such as quantum chemistry [l|, Q , 
field-linked states @ , the possibility for quantum compu- 
tation [1, [B[ , and long-range topological order Q . Fur- 
thermore, the closely spaced, opposite-parity internal lev- 
els of molecules, e.g. O-doublet, rotational, and vibra- 
tional levels, present new possibilities for precision mea- 
surement of fundamental physics @, H, S EJ EH • 

For these reasons, there has been much effort towards 
developing techniques to produce cold polar molecules. 
Current molecular cooling techniques can be character- 
ized as either association of ultracold atoms or direct 
cooling of molecules. While the association of ultracold 
atoms, either via a Feshbach or optical resonance 
[HI, is capable of producing molecules near quantum 
degeneracy [HI, these methods restrict experiments to 
a limited class of molecules - namely, those composed 
of laser-cooled atoms. Conversely, direct cooling tech- 
niques, such as buffer gas cooling [l5[ and Stark decel- 
eration fl6| are capable of producing cold samples from 
a wide range of polar molecular species. However, de- 
spite much work, the molecular density currently attain- 
able via these methods is limited to < 10 7-8 cm~ 3 . In 
the case of buffer gas cooling, the molecular density is 
limited by the presence of the helium buffer gas, which 
prevents further cooling; for Stark deceleration the poor 
efficiency (~0.01%) of decelerating molecules from a su- 
personic beam is the limiting factor [l7j . Furthermore, 
the 'single-shot' nature of these techniques makes it dif- 
ficult to accumulate cold molecules over time. 

Here we propose a new technique for producing contin- 
uous beams of cold polar molecules. When coupled with 
irreversible trap- loading [H|,[ll| this continuous beam de- 
celeration technique circumvents the poor efficiency as- 
sociated with traditional beam deceleration by allowing 




'Electronic address: eric.hudson@ucla.edu 



FIG. 1: (color online) Schematic of continuous beam decel- 
erator. Weak-field seeking polar molecules exiting the beam 
source are decelerated as they move into the high electric field 
region between the electrodes of the deceleration stage. At 
the stage center, STIRAP is used to transfer the molecules 
to either a less polarized weak-field seeking state or a strong- 
field seeking state (inset), resulting in net deceleration of the 
molecular beam. After the deceleration stage, the molecules 
are returned to their initial internal state via the inverse STI- 
RAP and the process repeated. For clarity, we have omitted 
the low-voltage electrodes between the deceleration stages, 
which are used to maintain molecular orientation. 



the accumulation of molecules over many seconds, simi- 
lar to the loading of a magneto-optical trap (MOT) [20| . 
Moreover, this method appears to allow the accumula- 
tion of densities appropriate for evaporative cooling and 
may represent a viable route to quantum degeneracy for 
a wide-range of molecular species. In the remainder of 
this work the method for decelerating continuous molec- 
ular beams is detailed and supported with calculations 
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of the overall efficiency. 

Because continuous beams are extended in space they 
are not amenable to traditional Stark decelerators, which 
rely on well-timed electric field pulses and a spatially 
compact pulsed molecular beam [16( . Figure [T] shows a 
method of deceleration that instead uses static electric 
fields, coupled with appropriately placed transition re- 
gions, to control the molecular internal state. A molecule 
in a weak-field seeking state (dipole moment anti-aligned 
with the electric field) exiting the beam source will be 
decelerated as it moves into the high electric field region 
between the electrodes. If at the center of this deceler- 
ation stage the molecular internal state is changed to a 
less polarized weak-field seeking state, the molecule will 
exit the deceleration stage with less kinetic energy than 
when it entered. After exiting the stage, the molecule 
internal state can be returned to its original state and 
the deceleration process repeated. These transitions can 
be driven with near unit efficiency by the technique of 
stimulated Raman adiabatic passage (STIRAP) }2lj . 

This approach has several advantages over traditional 
Stark deceleration. First, this technique allows for con- 
tinuous guidance of the transverse molecular motion dur- 
ing the deceleration process; transverse motion has re- 
cently been shown to be the limiting factor in traditional 
decelerator efficiency [22| . Second, as shown in the in- 
set of Fig. [TJ if the molecule's internal state is instead 
changed to a strong-field seeking state (dipole moment 
aligned with the electric field) at the electrodes center, 
the molecule will be decelerated as it moves out of the 
high electric field region. Since strong-field seeking states 
typically experience the largest Stark shift, this technique 
makes it possible to decelerate a beam of molecules to rest 
with substantially fewer deceleration stages than required 
in a traditional decelerator. Third, because molecules can 
be decelerated in strong-field seeking states and trans- 
versely refocussed in weak-field seeking states, this tech- 
nique can be more efficient for the deceleration of certain 
heavy molecules than the alternate-gradient deceleration 
technique 23]. Fourth, because this technique can de- 
celerate continuous molecular beams it is amenable to 
the use of a buffer gas beam source [24| , which provides 
a higher molecular flux at a much lower initial beam 
speed than room temperature molecular beams; this is 
extremely advantageous since the required number of de- 
celeration stages scales with the square of the molecular 
beam velocity. Finally, because this technique produces a 
continuous beam of decelerated molecules it is well-suited 
for dissipative trap loading, either into a static trap (l8l | 
or the recently proposed molecular MOT (l9j . 

For definitcness, consider a molecule with 1 E + ground 
and excited states, e.g. any compound containing one al- 
kaline earth atom and one chalcogen atom. The Hamil- 
tonian for the lowest vibrational level of each electronic 
state in the presence of an electric field is: 

H = Ti ;V + Bi^J 2 + p, ■ Edc, (1) 
Ti tV is the vibronic term energy operator with eigenfunc- 
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FIG. 2: Stark shift of 1 £+ states, labeled by k, and STIRAP 
transfer scheme. Molecules enter the deceleration stage in 
the k = 2 state and exit in k — 0. If the STIRAP transfer 
takes place at fixE/Bi = 15 then ~3.1 cm -1 is removed per 
deceleration stage. The figure inset shows the relevant Honl- 
London factors for SrO as a function of electric field. 



tions \ipi,v) where i — X,A,..., v — 0,1,2..., and Tx,o = 0; 
Bi tV is the rigid-rotor rotational constant for the ith elec- 
tronic level with vibrational excitation v; J is the to- 
tal angular momentum of the molecule; fl is the electric 
dipole moment operator; and Eoc is the externally ap- 
plied electric field. The eigenvectors of this Hamiltonian 
are \^i }V ,k) — J2j a j fc \J> M) \ipi, v ) with energies as shown 
in Fig. [2] as a function of electric field. Here the eigen- 
vector label k is equal to the value of J that describes 
the zero-field wavefunction. 

For the deceleration scheme considered here, the 
|*A,i,i) *~ l*x,o,2> and |*a,i,i) -> |*x,o,o) transitions 
are used to adiabatically transfer population during the 
deceleration process. If deceleration by exclusively weak- 
field seeking molecules is desired (for transverse guid- 
ing purposes) the |^x,o,o) state can be replaced by the 
|^x,o,i) state. Note that a two-photon transition be- 
tween states with |AJ| = 1 is possible here because the 
static electric field mixes states of opposite parity. In the 
rotating frame, the Schrodinger equation for the evolu- 
tion of the population of the relevant molecular states is 
given as [211 ]: 
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where T is the natural linewidth of the A state in radians 
per second, O n (i) = d ""^^ t - ) with E n (t) the envelope of 
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the laser's electric field amplitude, and c* k is the prob- 
ability of the molecule being found in the state |^i )t)) fc). 
The transition dipole moments are given in terms of the 
dipole operator, d, as 



d„ = (* A ,l,ll d |*X, 



0,2(2- 



-n)) , 



(3) 



which can be calculated from the relevant tran- 
sition moments: (^A,v| (J + 1, M| d |JM) \tpx,v) = 
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tains both the electronic and vibrational contributions 
to the transition dipole moment. 

As a specific example, Eqns. [1] and [2] have been 
solved for the strontium monoxide molecule, SrO, (By = 
0.33798 cm" 1 fi x = 8.9 D, B A = 0.3047 cm" 1 [25[, 
ua « 2.5 D H, T = 2ttx3.7MHz [27], and df$ = 0.43 D 
IH) in a 30 kV/cm electric field (n x E/B x « 15). For a 
geometry similar to Ref. [29[ , with laser beams of axial 
and transverse e -2 intensity waist wa x w_l of 1 x 10 mm 2 
and total power of 1 W, the transfer efficiency has been 
calculated to be > 90% throughout the deceleration pro- 
cess - because the pulse area scales as w" 1 , the transfer 
efficiency increases with deceleration stage number, N, 
as A 1 / 2 . At this field, and with transfer from the |2,0) 
to |0,0) state at the stage center, ^3.1 cm" 1 is removed 
from the molecular kinetic energy at each deceleration 
stage, allowing a 100 m/s beam of SrO to be decelerated 
to rest with only twenty deceleration stages. 

If the decelerated molecules are loaded into the mi- 
crowave trap of Ref. [3(| via a spontaneous emission load- 
ing scheme similar to that of Ref. [IH, then molecules 
exiting the decelerator with a final velocity less than the 
maximum trappable velocity, Wtrap, can be accumulated 
in the trap. Since, all molecules lose the same amount 
of energy AE per deceleration stage these trappable 
molecules enter the decelerator with a velocity in the 
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gas beam source can be described by the longitudinal 
velocity distribution, 
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the fraction of trapped molecules is then given as 
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where rjs.E. is the efficiency of the spontaneous emis- 
sion trap loading [l8[ and rj± is the fraction of molecules 
that remain within the transverse bounds of the decel- 
erator during the deceleration process [17j . For loading 
of the microwave trap by excitation of the I^a.i.o) <~ 
l^x.0.2) transition followed by spontaneous emission on 
the I^a.I.o) - ► |*x,o,o) transition [3l|, and assuming ex- 
citation saturation, Eq. [3] yields tjs.e. = 0.4 32]. For 
reference, Eq. [5] is plotted in Fig. [3] for SrO (vq = 100 
m/s and Av = 25 m/s) and the microwave trap param- 
eters of Ref. [3(| (wtrap = 2 cm -1 ). 

The transverse efficiency t]± depends sensitively on the 
details of the decelerator design and molecular beam, and 
can only be estimated from detailed Monte Carlo trajec- 
tory simulations. Since most of the deceleration takes 
place in a strong field seeking state and the strongest 
fields are near the electrode surface, molecules are defo- 
cused transversely as they are decelerated. Thus, a mech- 
anism for transverse refocusing of the molecular beam 
must be added to the decelerator to prevent complete 
loss of the beam due to transverse defocussing. The 
techniques of transversely focusing molecular beams are 
the subject of much on-going work [22l . |33| and a com- 
plete optimization of the bcamline is beyond the scope of 
the present work. Nonetheless, Monte Carlo trajectory 
simulations were performed for a decelerator geometry 
with a hexapole focusing element [34[ (t = 1 cm and 
V = 35 kV) added after every fifth deceleration stage. 
The results of these simulations are shown in Fig. [3] 
as open circles for a beam of SrO molecules with vq — 
100 m/s, Ad = 25 m/s, and transverse beam tempera- 
ture Xj_ = 4 K. At twenty stages of deceleration, 77 reaches 
a maximum of ~0.002. Thus, with an initial molecular 
beam flux of 10 13 molecules/s [24| and trap volume of 
V ~ |7r cm~ 3 30], the trap density accumulates at a 
rate of K w 5 x 10 9 cm~ 3 s _1 . For a background gas 
limited lifetime of r w 100 s, an average trap density of 
p = Kt w 5 x 10 11 cm~ 3 is expected. 

A recent proposal has shown how certain molecules 
may be laser-cooled [l!| and collected in a MOT. While 
most molecules are not amenable to the molecular MOT 
cooling scheme, many molecules possess relatively diag- 
onal Franck-Condon factors allowing the possibility of 
scattering many photons before optical pumping to un- 
desired states occurs. Thus, cooling only the transverse 
degree of freedom of a molecular beam, which requires 
the scattering of a maximum of a few thousand photons, 
may be possible for many molecular species. For refer- 
ence, simulations were also performed for a SrO molecu- 
lar beam that has undergone transverse laser cooling to 
T± = 10 fiK. As shown in Fig. [3] as closed circles, the 
reduced transverse beam temperature makes it possible 
to recover the maximum trapped fraction as predicted by 
Eq. [5] Of course, it may be possible to achieve the max- 
imum efficiency, even for Tj_ = 4 K, with more elaborate 

; transverse guiding schemes (e.g. addition of more focus- 
ing elements, different order focusing elements, voltage 

(5)scaling 22]) and this is the subject of on-going work. 
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FIG. 3: Results of three-dimensional Monte Carlo simulation 
of the deceleration process. The fraction of molecules that 
exit the decelerator and can be coupled into the trap is dis- 
played versus the number of deceleration stages used. Open 
circles are for a traditional buffer gas beam of SrO molecules 
(Tj_ = 4 K). Closed circles are for a molecular beam that has 
undergone transverse laser cooling (T± = 10 /dK). The line 
represents the maximum trappable fraction as predicted by 
Eq. E 

An isolated dilute gas can be efficiently cooled by self- 
accelerating evaporative cooling as long as the product of 
the elastic collision rate and trap lifetime exceeds ~150 
[35|. For polar molecules in their absolute ground state, 
the collision rate is expected to be accurately predicted 
by the Eikonal approximation for T > 200 fK [13, GsB] 
and is given as [13] as 7 = pcrv = w 2 x 10 5 Hz. 
This large scattering rate (by atomic physics standards) 
is due to the long-range nature of the dipolar interac- 
tion between polar molecules in their absolute ground 
state. Since jt w 10 7 , the criteria for efficient evapora- 
tion is clearly satisfied and evaporative cooling appears 
to be a viable route to quantum degeneracy. However, 
the anisotropic dipolar interaction may lead to a signifi- 
cantly peaked differential cross-section, which, due to the 
reduced momentum transfer for small angle scattering, 
could reduce the efficiency of the cooling below this esti- 
mate. Furthermore, species dependent, density limiting 
effects such as three-body recombination have not been 
studied in detail and could lead to unexpected losses; in 
analogy to evaporative cooling of cesium it is likely that 
an externally applied electric or magnetic field can be 
used to mitigate the effect. 

Because the production of the necessary electric fields 
and ultra-high vacuum is now routine, the main technical 



challenge in constructing the continuous beam decelera- 
tor is the required voltage stability. Since the two-photon 
detuning, A, depends on the electric field, any change 
in the field during the deceleration process diminishes 
the efficiency of the STIRAP process. At 30 kV/cm 
the energy dependence of the 1 20) \ipx,v) state on elec- 
tric field is weak compared to the field dependence of the 
1 00) |V>X,v) state. Therefore, A depends on electric field as 
A w (no) S\E\, where (no) is the ground rotational state, 
laboratory-frame dipole moment (6.9 D for SrO (30| ) and 
8\E\ is any drift of the electric field during the decelera- 
tion process. Maintaining maximum STIRAP efficiency 
requires A < 2itv/wa- Since the STIRAP efficiency is 
oc ^/wa/v, while the pulse time is sa wa/v, the beam 
waist may be decreased with the square of the velocity to 
maintain maximum transfer efficiency. Thus, the tightest 
constraint on field stability comes from the first deceler- 
ation stage, which in the configuration of Fig. [2] leads to 
A < 2-7TX100 kHz, or a required electric field stability of 
<5|-E|/|i?| ~ 10~ 6 . For instance, a 1 cm electrode spacing 
requires a challenging, but possible, voltage stability of 
0.1 V out of 30 kV. Any inhomogeneity between deceler- 
ation stages may be offset by local voltage adjustments. 
Moreover, it may be possible to utilize a field inhomo- 
geneity to produce an adiabatic rapid passage transfer 
|38| , greatly reducing the constraints on field stability. 

In summary, a general technique for producing cold po- 
lar molecules has been described. This method promises 
large improvement over other direct cooling methods by 
allowing deceleration of a continuous molecular beam as 
opposed to the single-shot nature of current techniques. 
By coupling the decelerated molecular beam into a trap, 
densities > 10 11 cm~ 3 seem attainable, which represent 
an improvement of up to 4 orders of magnitude over cur- 
rent technology. At these densities it appears efficient 
evaporative cooling to quantum degeneracy may be pos- 
sible. Finally, in this work STIRAP has been used to 
accomplish the necessary state transfer, however, in the 
presence of an electric field, parity is mixed and a single 
photon can transfer population between the relevant lev- 
els. Thus, it may be possible to perform direct transfer 
by THz radiation, further simplifying the design. 
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